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1 Introduction 



Doping of antiferromagnetic Mott insulators causes frustration which has a profound effect on 
the magnetic properties of such systems. High temperature superconductivity, charge ordering 
and anomalous transport properties have been observed in doped transition metal oxides with 
perovskite structure (see Ref. for a review). Numerical studies of one-dimensional models 
proposed for doped Nickel-oxides and Manganites for such systems showed strong tendencies 
toward ferromagnetism and phase separation. For a better understanding of these phenomena 
requires to take into account strong electronic correlations. A commonly used starting point 
for a description of such systems is the ferromagnetic Kondo Hamiltonian 

n = -tJ2 {4aCja + h.c}j - J^f ^ (Ti ■ S, (l.l) 

(ij)a i 

where the first term is the kinetic energy given in terms of canonical fermionic creation and 
annihilation operators cj^ and Cjo- for the itinerant electrons on site i with spin a =t, | and the 
second is the ferromagnetic Hund's rule coupling between the spins af = J2ai3^la ('^'^)a/3 ^*/3' 
k = x,y, z, of the itinerant electrons to a localized spin S„. A Coulombic repulsion to suppress 
double occupancy in the itinerant band is implicit. Large Hund coupling Jh favours the align- 
ment of the itinerant and the localized spin, i.e. spin eigenstates with the maximum allowed 
total spin. Hence, for local spins of length S — 1/2 the degrees of freedom that need to be kept 
for an effective theory of the low lying modes in this model are spin-S* "spins" and spin-S* — 1/2 
"holes". Within the double-exchange Hamiltonian classical "background spins" Sj are 
used to approximate these holes. However, the non trivial phases arising in a full quantum 
mechanical treatment of these spin degrees of freedom may induce low-energy modes which are 
essential for an understanding of the magnetic properties of these systems [Q. 

To obtain an effective lattice model on the 45* + 1 dimensional local Hilbertspace one has to 
eliminate the other allowed spin configurations in a perturbative analysis for Jh ^ t 0,0. To 
leading order in Jh the effective Hamiltonian of the resulting model is simply the projection of 
(|1.1|) onto the states listed above. 5'f/(2)-invariance implies that this operator can be written 
as a polynomial of spin-operators: 

Kff ~ -tJ2v,,Qs{yij) = S, ■ S,/S{S - 1/2) (1.2) 

where the operator Vij permutes the states on sites i and j thereby allowing the holes to 
move. For large but finite Jh one finds additional antiferromagnetic Heisenberg exchange terms. 
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Hamiltonian operators of this type have been used as a starting point for studies the phase 
diagram of doped transition metal oxides by numerical diagonalization of small clusters 

In this paper we introduce a class of exactly solvable models in one spatial dimension which 



generalize the supersymmetric t-J model and a model for doped spin-1 chains |T0|. In 
spite of the appearance of several additional couplings guaranteeing integrability these models 
may provide further insights into the peculiar properties of these compounds. In the following 
section we shall use the framework of the Quantum Inverse Scattering Method (QISM) to 
construct families of vertex models making use of so called 'atypical' representations of the 
super Lie algebra gl{2\l). The spectra of the corresponding commuting transfer matrices are 
obtained by means of the algebraic Bethe Ansatz (BA). In Sect. |]we derive local Hamiltonians 
(i.e. operators involving nearest neighbour interactions on the lattice only) similar to the ones 
discussed above using a fusion procedure. In Sec. |^ integral equations determinining the spectra 
and thermodynamic properties of these models are obtained from the BA equations. From these 
equations we obtain the phase diagram of the doped spin chains in a magnetic field for low 
temperatures T <^ H and the low temperature properties for vanishing field = in Sects. ^ 
and ^ We conclude with some remarks on a possible 5'[/(2)-invariant effective field theory 
description for the low energy/low temperature sector of these systems. 

2 Construction of the models 

Below we will construct a class of integrable vertex models from solutions to the Yang Baxter 



equation (YBE) which are invariant under the action of the graded Lie algebra gl(2\\) |TT, 
The nine generators of (7/(211) are classified into even (1, 5*^, S*^, B) and odd (V±, W±) ones 
depending on their parity w.r.t. grading. The even generators are the spin operators S*" form a 
SU{2) subalgebra with commutation relations [S*"^, S'^] = ±5*^, [5*+, 5^] = 2S^ and the U{1) 
charge B which commutes with the 3°": [5, S^] = [B, S^] = 0. The commutators between 
even and odd generators of the algebra are 

[S\V±] = ±iv±, [S^,V±] = 0, [S^,V^] = V^, [B,V^] = 

[S\W±] = ±^W±, [S^,W±] = 0, [S^,W±] = W^, [B,W±] = -^W±, (2.1) 

and the odd generators satisfy antzcommutation rules 

{V±, V^} = {V±, V^} = {W±, W^} = {W±, W^} = 0, 

{V^±,W^±} = ±V ' {V±,W^} = \{-S^±B). (2.2) 



The 'typical' representations [b, s] of this algebra can be characterized by the eigenvalues of 
operators B and on the multiplet with largest total SU{2)-spm [|Tl|,[12|. Their dimension is 
85" and they can be decomposed into two spin-(S' — 1/2) multiplets with charge 6 ± 1/2 and 
a spin-S" and a spin-(S' — 1) multiplet with charge b with respect to the S'f/(2)-subalgebra of 
gl{2\l). In the following we shall be particularly interested in the (45* + 1) -dimensional so- 
called 'atypical' representations [5*]+ which contain two multiplets of spin 5* and {S — 1/2) and 
corresponding charges b = S and b = S + 1/2. Choosing a basis {\b, s, m)} in which B, and 
are diagonal, the nonvanishing matrix elements of the fermionic operators are 



{S + -^,m±^\V±\S,S,m) 
{S,S,m\W±\S+^,S-^,mTl) 



± 



S ^ m 



S ±m 



(2.3) 



Tensor products of atypical representation can be decomposed as 



[sr = [s + s'i 



S + S' + ^,S + S'-^ 



S + S' + ^,S + S'~^ 



s + s' + ^,\s-s'\ + ^ 



(2.4) 



The irreducible components in this tensor product can be identified by the action of the 
quadratic Casimir of the algebra 



B' - W-V+ + W+V- - V-W+ + V+W. 



(2.5) 



which has eigenvalues on and (s^ — 6^) on [6, s\. 

Choosing an irreducible d-dimensional representation oi gl{2\l) acting on a quantum space 
V ~ C^, it is straightforward to verify that the ^-operator 



/ ii + 2iB 



£(/i) 



iy/2W. 



iV2V+ fi + i{B + S') -iS+ 
\-iV2V- -iS- fi + i{B-S')J 



(2.6) 



written as a matrix in the three-dimensional matrix space Ai solves the Yang-Baxter equation 



(2.7) 



with the i?-matrix 



R{\) =b{X)I + a{\)U , a(A) 



A 



6(A) 



A + i 



(2.8) 
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Here / and 11 are the unit and graded permutation operator acting on the tensor product 
Ail ® A^2 of matrix spaces in which /^-operators act in ( |2.7] ). Assigning Grassmann parities 
Cj G {0, 1} to the basis of these spaces the matrix elements of 11 are 

= i-^r''''SnrA.n ■ (2.9) 

Similarly, the matrix elements of the operators acting on the tensor product of these spaces pick 
up signs (A (g) 5)*^^.^ = (—ly^ii'^^i+^n^Aij^jj^Bi^j^ due to the grading of the basis. Considering 
the /^-operator as a linear operator acting on the tensor-product of spaces Ai ® V with the 
fundamental three-dimensional representation [1/2]^ acting on the matrix space its gl{2\l)- 
invariance can be established by rewriting (|2.6|) as /i — 1X2 in terms of the Casimir operator 
(|2.5|) on the tensor product (up to a shift of the spectral parameter). 

The intertwining relation ( |2.7| ) implies that the monodromy matrix, defined as the matrix 
product 

T(A) = £i(A)/:i„i(A)---£i(A) (2.10) 



of £„-operators ( ^.61 ) with entries acting on different quantum spaces V„ satisfies a Yang-Baxter 
equation with the same i?-matrix (|2.8|): 



R{\ - /i) (T(A) ® T(/i)) = (T(/i) ® T(A)) i?(A - /i) . (2.11) 

As an immediate consequence of this identity the transfer matrix given by the matrix super 
trace of T 

commutes for different values of the spectral parameter fi thus being the generating functional 
for a family of commuting operators on the graded tensor product of L quantum spaces which 
we will identify below with the Hilbert space of an integrable spin chain. The subscript to 
the transfer matrix is used to label the dimension of the matrix space of the corresponding 
monodromy matrix. 

The spectrum of this transfer matrix is obtained by means of the algebraic Bethe Ansatz 
(ABA) Ijl^. As a consequence of the grading different sets of Bethe Ansatz equations (BAE) 



follow from different orderings of the basis [13]. We now restrict ourselves to representations 
[S]^ in the quantum spaces V„ where we choose the state |0)„ = 15, 5*, S)n as our reference state. 
The two other possible sets of BAE for this model are given in Appendix ^, their equivalence 
is shown in Appendix ^ 
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Since the ABA for the transfer matrix ( p.l2|) is completely analogeous to the case considered 



in |T^,|1^,0, we only sketch the main steps leading to the BAE and the spectrum. The action 
of (U-6\) on the reference state is 



/ iJ + 2iS 







\ 



|0), 



(2.13) 



fi + 2iS 
\ -tV2V- -iS- /i ) 

Similarly, acting with the monodrony matrix ( p.lO| ) on the state \VLs) = [0)^ (S) • • ■ ® |0)i we get 

/ {fi + 2iS)^ \ 

Tifi)\ns)= (/i + 2z^)^ l^s) . (2.14) 

\ Ci(/i) C2(/X) /i^ / 

Hence, \^s) is an eigenstate of the transfer matrix ( ^.12| ) with eigenvalue {—fi^) (we have chosen 
the grading ei = 0, e2 = €3 = 1 in the matrix space here). The operators Ci(A) and C2(A) 
create a hole and lower the spin of the system respectively. For eigenstates of ^35 (A) with Nh 
holes (generating sites with spin 5* — 1/2) and magnetization = LS — \Nh — we make 
the Ansatz 

|Ai, . . . , A„|F) = C„,(Ai) ■ ■■CaX^n)\ns)F"^ (2.15) 

where n = Nh + Ni. Using the algebra of the operators in ( p.llD we are led to an eigenvalue 
problem for the amplitudes F"'" '"'^ which is solved by a second Bethe Ansatz parametrized by 
'hole rapidities' {i^a}a=i- Finally, we find that ( |2.15|) is an eigenstate of ( p.l2| ) with eigenvalue 



fi-Xj 



+ (/i + 2iS)'' Yl 



a=l 



1- n 



Aj — fi + t 
Xj — fj, 



(2.16) 



provided the spectral parameters Xj = Xj — iS and Ua = — iS + i/2 satisfy the following set 
of BAE 



Xj + iS 
Xj — iS 



n 



Xj - Xk + i 
Xj — Xk — i 



n 



^Xj - Va + ^ 



J 



,Nh + N 



(2.17) 



Nh+N 

n 

k=l 



l^a — Ak + 
— Xh — 



a = l Nh . 
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3 Doped spin chains 

Choosing the fundamental three dimensional representation [1/2]^ of gl{2\l) in ( |2.6|) , the C- 
operator taken at fi = —i becomes a graded permutation operator on the tensor product Ai^V 
of matrix and quantum space. Hence the transfer matrix ( p.l2[ ) generates a translation by one 
site on the lattice for this value of the spectral parameter. Expanding the logarithm of tslfi) 
about this shift point we can construct a Hamiltonian with nearest neighbour interactions only 

= ±... (3.1) 



■ 9 , , , 
-.^lnt3(/.) 



M=-« n=l 

which is the supersymmetric t-J model (see e.g. [1^). In this case, Eqs. ( |2.17|) are known as 
Sutherland's form of the BAE for this model pj. 

For representation different from [1/2]^ it is not possible to construct a local Hamiltonian 
directly from the ^-operators ( p.6| ) (they can be used to construct t-J models perturbed by 
integrable impurities though [p!7|-[T^). To obtain a homogenous lattice model such as ( |3.1| ) new 
^-operators have to be found which act on tensor products of matrix and quantum spaces of 
the same dimension. Their (7/(2|l)-invariance implies that they can be expressed as sum over 
multiples of the projectors on irreducible components of the tensor product of representations 
in the two spaces. Noting that the spin multiplets at charge {S + S') in the tensor product ( p.4| ) 
are just the ones obtained by adding two spins of length S, S' (and similarly spins {S — 1/2), 
(5" — 1/2) at charge (5 + 5" + 1)) we find linear operators C^^^'^dx) acting on spaces carrying 



atypical representations [S]_^ and [S']_^ from the 5f/(2)-invariant ones constructed in Ref. ||20 
namely: 

S+S' ., S+S'-l m _ 

fc=|S-5'|+l ^ m=\S-S'\ k=\S-S'\+l ^ 

Here Va is a projector on the (7/(21 l)-multiplet A in the tensor product [5']_^ [S']^. Choosing 
one of the representations to be [1/2]^ and comparing this expression with (|2.6| ) we find 

£{1/2,5} ^ - (^ - ^{S + 1/2)) . (3.3) 

fi + i[b + 1/2) 

The new /^-operators satisfy the YBEs 

i?5.5.(A - fi) (C^'^'^H^) ® C^'^'-^Hfi)) = (C^'^'^Hf^) ® C^'^'^HX)) Rs,s,{\ - ^) (3.4) 

where Rss' = IIC^^^'\ As a consequence of this relation the transfer matrices of the corre- 
sponding vertex models 

t^^°^}(/x) = sTro [C['"'\fi) ■ ■ ■ C\'^''\f^)) (3.5) 



(the product of /^-operators and the super trace are taken in the (45*0 + l)-dimensional matrix 
space) commute, i.e. [t'f'^°'^^(/i), f'^'^^'^^(A)] = for all So, Si. From ( |3l^ ) we observe that 
C^^{fi = 0) oc n. As in (|3.1| ) an integrable Hamiltonian with nearest neighbour interactions on 
the lattice can be constructed by taking the logarithmic derivative of \iit^^^^{n) at /i = 0. 

The eigenstates of the transfer matrices t^^'^^{iJ,) are parametrized by the roots of the BAE 
(|2.17|) . To compute their eigenvalues we need the so-called fusion relations between these 
operators for different S' which are obtained from considering tensor products of different 
matrix spaces. 

Starting with the YBE ( |3.4|) for 5*1 = 5*2 = 1/2 we observe that choosing X — fi = —i the 
matrix Ri i is proportional to a projector onto the five- dimensional subrepresentation [1]^ of 
the tensor product [1/2]^® [1/2]^. This implies that the ^-operator £(/i) = C^^/"^'^^ (^ — 
£{1/2,5} -|_ 2/2) satisfies the condition 

P[l]^£(/i)P[3/2,l/2] = . (3.6) 

Consequently, it can be rewritten as 

C{li) ~ (3.7) 

1^ £11-^1 (^) ) 

after a proper reordering of the basis in Adi ® M.2- Here £'f[^/^'^/^]''^J'(yu) is a 4 x 4 matrix 
acting on the [3/2, 1/2] component of this tensor product. Building a monodromy matrix from 
L copies of (|3.7|) we obtain the fusion relation for the corresponding transfer matrices 



t(/i) = sTr [CLi^i] ■ ■ ■ 

= t^y^^s} _ ^/2) ^{1/2,5} ^ ^/2) = t^l'^>(/i) + t{[3/2,l/2],5}(^^ (3_8) 

and an equivalent equation for their eigenvalues A^'^^{fi). Since there are still two unknown 
functions of fi on the RHS of this equation it is not possible to determine the spectrum of the 
new transfer matrices from ( p.8|) alone. As additional information we make use of the fact that 
the eigenvalues of the transfer matrix are analytical functions of fi, i.e. the residues at their 
simple poles vanish as a consequence of the BAE ( |2.17| ). Complemented by the trivial action 
of t^^'^^fi) and t{[3/2,i/2],5}(^^^ ^j^g pseudo vacuum 1^5) this allows to compute A'"^^''^^ with 
the result (see e.g. [pl[] ) 



I fi - tb fi + t{b - 1)\ T-r - Aj + 2 TT- - Z/„ - 2 K yr Xj - fi + 2t 

\fx + tS fi + tiS+l)J J-jL /i - A,- A A - zy, - f I " n A,-/i + 2 

^^f^§ n n ^ (3.^) 

-t- il^O i;/ Aj ^ Q,=l '^a 2 I j=l J ^ ^ I 

As observed above, the local vertex operator C'^^'^^^fj, = 0) becomes a graded permutation 
operator on the tensor product of the five-dimensional matrix space and the quantum space in 
which the representation [S = 1]^ is acting. Hence we can proceed as for the the case of the 
fundamental representation above and obtain the local lattice Hamiltonian for a spin-1 chain 



doped with 5 = 1/2 holes introduced in Ref. 



n^'^ = -^^\n{ti'^'H^^)) - SL = {^^'^-,1 + Kjli] " • (3-10) 

M=0 n=l 

The exchange and kinetic part of the Hamiltonian expressed in terms of spin operators Sj with 
Sf = Si{S, + 1) with Si = lor 1/2 read 

= I (^^S. ■ S, - 1 + 6s^s,,^ (1 - (S. ■ S,)^)) , 

-^^r = - (1 - ^s.,s,) (S. ■ S,) . 



Vij permutes the spins on sites i and j. The corresponding eigenvalues of ( 3.10 ) are obtained 
from ( p.9|) : adding an external magnetic field H and a chemical potential 

(we have added an external magnetic H field and a (hole) chemical potential fi to the Hamil- 
tonian) . 

To proceed to higher S we iterate the procedure used above: for 5*1 = 1/2 and S2 = S' — 1/2 
arbitrary we use the fact that R^ s'-^if^ = —iS') oc "^[5']+ in the YBE (|3.4| ). This leads to the 
fusion relation 

which allows to determine the eigenvalues of t^^''^^{fi) from the known ones of t^^^"^'^^ (fi) and 
^{S' -1/2,3} as 



S+S' \ ^ N^+Ni 



A'-^'(^|{^'*.'^ - n ^ n 

\k — \S .51+1 / j 1 



+ [^i-^{S-S')y{...} . (3.13) 



The terms in braces are determined by the fusion equations together with the vanishing of the 
residues at the simple poles of A^^''^^{fi) due to the BAE ( |2.17| ). They do not contribute to 



the spectrum of the nearest neighbour spin chain Hamiltonian 

d 

n'^^^ = -z — In [t^^^^\^i)) + const. (3.14) 

whose eigenvalues are 



5: 

k=\ ^ K / ^—i \ / 



4 Thermodynamic Bethe Ansatz 



To study the thermodynamics of the doped spin chains ( 3.14| ) we have to analyze the BAE 



( |2.17| ). In the thermodynamic limit L —>■ oo general solutions are known to consist of real hole 
rapidities Ua and complex n- strings of spin-rapidities 

Xf = X(p^+' (n + l-2k) , k = l,...,n. (4.1) 



Now we consider solutions of ( |2.17|) built from hole rapidities and M„ A-strings of length 



n. Rewriting the BAE in terms of the real variables and and taking the logarithm we 
obtain 

(^«) = ^ , zi-^ [Xf) = , (4.2) 
where jj"^ and are integers (or half-odd integers) and the functions Zi are given as 

^ oo Mm Nh 

2nzi-\X) = e„,25(A)--$^$^S_(A-Af))+-5^e„(A-z.„) 

m=l j=l a=l 



oo AIn 

2TTZr(iy) 

ij 

n=l j=l 



^ oo AIn 

rEEs-f-^r) («) 



where 0„(x) = 2 arctan(2x/r;,) and 

0nm{x) = ^m+n-1 i^) + i^) + . . . + 0\m~n\+l i^) y 

(x) + . . . + 2^1 (x) + (1 (x) . (4.4) 
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The quantum numbers jj"'' and la in (^4.2|) uniquely determine a particular eigenstate of the 
system. The asymptotic behaviour of the functions ( |4.3| ) determine their possible values. This 
allows to introduce densities piy) of the hole rapidities, cr„,(A) of the A-strings and the corre- 
sponding hole densities p(z/) and (t„(A) with 

d d 
an{x) + cr„(x) = -Q^^^r\^) > P{^) + P(a;) = -q^^c{x) . (4.5) 

In the thermodynamic limit L — oo with N^/L and M„/L held fixed these equations become 
linear integral equations 



an{x) = {An,2S * s) (x) - ^ {A^m * Oni) {x) + (a„ * p) [x] , 

m 

p(x) + p(x) = ^ (a„ * (T„) (x) . (4.6) 

n 

Here, (/ * (?) [x) denotes a convolution, 27ra„(x) = Q'Jyx) = 4n/(4x^ H-'^^), s{x) = 1/(2 cosh vrx) 
and 

Anmix) = — S^„(x) + SnmSix) . (4.7) 

Ztc 



Similarly, the energy (|3.15| ) in the thermodynamic limit can be rewritten in terms of the densities 
E/L = f2j (ef (a;) + riH) (t„(x) - j dx (^fi + ^H^ p{x) (4.8) 

where en\x) = —271 {An,2S * s) {x) are the bare energies of the A-strings. 

At finite temperature the equilibrium state is obtained by minimization of the free energy 
F = E — TS by variation of cr„ and p. Here S is the combinatorical entropy ||22 | 

oo „ 

S/L = ^ / dx {(cr„ + a„) In (cr„ + cr„) - cr„lno-„ - o-„lncr„} 

n=l ^ 

+ y dx {(p + p) In (p + p) - plnp - plnp} . (4.9) 

As a result we obtain the thermodynamic Bethe ansatz (TEA) equations for the energies 

e„ = Tln(cr„/(T„) of A-strings and k, = Tln(p/p) for the hole rapidities 

T 
2^ 



^n{x) -T^Yl -nm * In [l + c"^™/^] (x) + Ta, * In [1 + e"'^/^] (x) = ei°) (x) + nH 



k{x) + T ^ a„ * In [1 + e"^'"/^] (x) = - U + ^i/ j (4.10) 
An alternative form of these equations can be obtained by using the identity J2k i^nk * Akm) (x) 

5nm 5{.x) with 

Cnmi^x) 6nm^(^x) ('^n+l,r?i + K~l,m) S{x) . (4.11) 
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This allows to rewrite the integral eqs. ( [4.6|) as 

{a2s*s){x) = p(x) + ([1 + 02]"^ *p) (x) + (s *(Ti) (x) . (4.12) 
Similarly, we find for the energy of the corresponding state 

E/L = E^Q^^/L- j dx[27r{a2s*s)+p]p{x)+ j dx27is{x)a2s{x) (4.13) 

— lim Hn / dxanix) 

( S) 

where Eq is the ground state energy of the spin-^* Takhtajan-Babujian chain in a vanishing 



magnetic field 



,(5) _ I - Efc=i 2fc3T for integer S 



-2 In 2 - J2tZi I for half-integer S 



Finally, an equivalent form of the TBA equations ( 4.1CI|) is 



e^{x) = T (s*ln [1 + e^"-^/^] [1 +e^"+^/^]) (s) 

-27r5„,25 six) - 5n,iT {s * In [1 + e"''/^] ) (x) , (4.15) 
subject to the condition lim„_,oo(en/''^) = H and 

- [27ia2s * s{x) + /i] - T (s * In [1 + e^^/^j ) (x) = k{x) + T (i? * In [l + e"''/^] ) (x) (4.16) 

where i? = 02 * (1 + 02)""'^. 

In terms of the solutions to these equations the free energy reads 
F/L = e'^^^/L -T j dxs{x) In [l + e^^^W/Tj _rj. j ^^^^^^ ^ ^-^^^^-^ ^ g-«(x)/Tj _ ^^^^-^^ 

5 Zero temperature phases in a magnetic field 

In the limit T the TBA eqs. ( |4.1CI| ) become linear integral equations. As a consequence 
of ( [4.15| ) only ei(x), €25 (x) and k(x) can take negative values for certain x. Hence we have 
to solve three coupled integral equations for these quantities which in turn determine all other 
dressed energies 

en(x) + ^ js^i * eS"^ + S;^25 * } (x) - a„ * ^^-\x) = (x) + nH 



K(XI 



ai*e[^ + a2s * } (x) = - ( p + ] , (5.1) 



where /^(x) = ^(±/(x))/(x). 

To discuss the solutions of these equations further we have to distinguish various cases: 
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5.1 ii>H/2 

In this regime we have k,{x) < and ei(x) < for all x. This allows to express these functions 
in terms of the remaining unknown function e2s{x). From { f).l\) we find 



K,{x) = — 2/i + * e[' 



(0) 



[X 



ei(x) = e!°^(x) - - - («25-i * 4s) (x) 

en>i{x) + ^ (s'„_i,25-i * 4's) (x) = -27r {An-i,2s-i * s) {x) + {n- l)H 

The last set of equations can be identified with the integral equations for the dressed energies of 
the spin-S* — 1/2 Takhtajan-Babujian model, hence this regime corresponds to the completely 
doped case (i.e. = L holes). For mag netic field H > //(^^Va) ^jth 

k=l 

the system is in a ferromagnetically saturated state with maximal magnetization = L{S — 
1/2). 

5.2 -H/2 < ii<H/2 

Here we find from (|5.1| ) that k(x) = k^^\x) < for all x while ei(x) can take non- negative 
values. Eliminating k,{x) from the integral equations for ei and €23 we obtain 

ei(x) + |a25-i * £25^1 (x) = ef'\x) - fi + ^H , 

f _ ^^'^ _ 1 45-1 
(^2s{x) + < 025-1 * eS"'' + 2 ^ a2A; * £25"* \ (^) = 4s (^) - ^ — H . 

In this regime we find ei > and €23 > for 

l^H-it-^,lH-l\ . (5.3) 

2 ^ 2A; - 1 ' 2 S f ^ ^ 

fc=i ) 

Positive dressed energies for the A-strings imply M„ = and from ( |4.6| ) we find that = 
in this region. Hence for T ^ ( |5.3| ) belongs to the ferromagnetically saturated phase of the 
undoped system, namely the spin-S* Takhtajan-Babujian model. This phase exists for magnetic 
fields H > H^^\ 

Increasing the hole chemical potential to values fi > H/2 — 2/S holes are added but the 
ground state continues to be fully polarized: For the dressed energies this corresponds to €25' > 
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H < mm 



and while the real spin rapidities fill all states with negative ei(x). As a consequence of the 
free fermionic nature of this state these dressed energies can be expressed in terms of their free 
values 

6i(x) = 6f)(x)-(^^-|) . (5.4) 

The lower boundary of this phase in the ^-H plane is given by the condition min^ {^2s{x)} = 0. 

For magnetic fields {AS — l/2)H < fi + 4 Ylt^i 1/(2^ — 1) the ground state is a filled sea of 
A-strings of length 25* with negative energy 

e2s{x) + I E * ey I (x) = eg(x) + (^2S - H - fx . (5.5) 

The other dressed energies k,{x) < and 6^^25(3;) > can be expressed in terms of the solution 
of this equation. In this region of parameters the system has a finite concentration of holes and 
overturned spins. Although one might naively expect two types of massless excitations in such 
situation only one such branch with dispersion (|5.5|) is realized in this system which turns out 



to describe the charge excitations. Hence spin excitations are gapped in this regime [24 



5.3 /i < -H/2 

Again we find several phases that can be characterized by the configurations of strings present 
in the ground state: For H > H^^^ all dressed energies are positive corresponding to completely 
polarized undoped state. 

For smaller magnetic fields e2s{x) takes negative values in some interval to be determined 
from 

e„(x) + ^ * 4^ } (^) = (^) + nH . (5.6) 

These are the Bethe ansatz equations of the spin-S* Takhtajan-Babujian chain. This phase 
becomes unstable against hole creation for chemical potentials 

> I da2six)ei-\x) -\h^4^ (^^) " ^ (^^) ^ = (5-7) 

{ip{x) is the digamma function). Beyond this line the ground state is built from a filled sea 
of A-strings with energies < and another sea of charge rapidities with energies k{x) < 0. 
Increasing the chemical potential further negative energy solutions for the real spin rapidities 
ei appear giving rise to a third condensate of Bethe rapidities. 

From the cases considered above we obtain the qualitative zero temperature phase diagram 
of the doped spin-S system in the fi-H plane presented in Fig. [l|(a). Using Eqs. ( [4.6| ) the 
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corresponding phase boundaries can be given as a function of the hole concentration x = Nh/L. 
For S = 1 this is shown in Fig. |l](b). 



6 Low-temperature thermodynamics at H = 

Further simphfication is possible in the case of finite doping in a vanishing magnetic field 
which corresponds to chemical potentials fi G [ip{{2S + l)/4) — ^^{{23 + 3)/4), 0]. Furthermore, 
ei(x) < 0, e2s{x) < for all x while k,{x) < for < Q in this regime. All other dressed 
energies vanish in this limit. Eliminating the e„ from the equation for the energy of the holes 
we obtain 

rQ 

- [2'Ka2s * s{x) + fi] = k{x) - dyR{x-y)K{y) (6.1) 

J-Q 

where the Fermi point Q is determined by the condition k{±Q) = 0. Similarly, the density of 
hole rapidities po{x) in this regime is given by 



rQ 

/ dyR{x - y)po{y) = 025 * s{x) . (6.2) 
J-Q 



Excitations with charge rapidities near ±Q are massless. The velocity of this charge mode can 
be obtained from the dispersion (|6.1|) 

1 Oko 



27rpo(Q) dx 



(6.3) 



Similarly, one has massless excitations near x = ±00 in the magnetic sector with energies ei {x) 
and €23 (x) with velocities 

f 25 = lim — -— = TT , t>i = lim -— = ^ . (6.4) 

x-*oo 2na2s{x) x^oo2nai{x) 2 J^^ dye'^y po{y) 

As a consequence of the behaviour of the dressed energies as — * we can replace k in 
Eq. ( [4.15| ) by its zero temperature value k,o{x) and the driving terms by their asymptotics to 
obtain the leading low temperature behaviour. As a result we get 



Ts * ln[l + e^"-i(^)/^] [1 + e^"+i(^)/^] - 2Tx5n,2S e""'"^' - 27rA(5„,ie-"l^l (6.5) 



where 2t[A = — J^^ dye^^ KQ{y) . 

To move further we have to separate the contributions to the free energy stemming from 
the charge-sector from those due to the e„. Considering low temperatures again the leading 
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contributions to k come from the vicinity of the Fermi wave vectors ±Q. In this region one can 
safely neglect contributions to Eq. ( |4.16| ) from ei and rewrite it as 

-[27Ta2s * s{x) +ij]-TR* ln[l + e'l'^^)'/^] 
rQ 

= - / dyR{x - y)n{y) (6.6) 



where Q is determined by the condition k{±.Q) = 0. Using the procedure introduced by 
Takahashi |2^, we can rewrite the free energy (|4.17|) as F/L = Eq^^ /L + fc + fs where 

fc = -T J dxpo{x)\n[l + e-^'''^''^^^^] ^ -nT^ev , (6.7) 

/, = -T J dxs{x)\n[l + e'^^^''^/^] -T j dx{s* po){x)\n[l + e''^='^/^] . (6.8) 

Now the thermodynamics is described by Eq. (6]7) representing a scalar bosonic mode (the 
charge sector) and by Eqs. (|6.5| ) and (|6.8|) for the spin sector. 

At low temperatures the spin contribution fs is dominated by contributions from the regions 
f25 exp(— 7r|x|) ~ T where |e25| ~ T and the second one by the regions Vi exp(— 7r|a;|) ~ T. The 
leading temperature dependence of fs at low T can be obtained by rewriting ( p.5[ ) for large 
positive X as 

M^) = s * ln[l + e^"-i] [1 + e^"+^] - Sn,2S e"'^" - A5„,ie--^ (6.9) 
in terms of the T-independent functions 

(Pn{x) = — e„, [x In 



T V 27r, 

In the low-T limit we can replace s{x) and s * po{x) in ( |6.8|) by their asymptotics to obtain the 
free energy 

/, (£Hi + £i) . (6.10) 



6 V^2S 

Such an expression is characteristic of a system two with massless excitations with velocities 
(|6.4|) . In cases where these excitations can be characterized by different observable quantum 
numbers the coefficients q are the central charges of the underlying Virasoro-algebra thus 
determining the universality class of the system. In this case they are given in terms of the 
solutions of (|6.9| ) by 



= ^ J da;e-""ln [l + e^^^^")] , ci = ^ J dx (Ae""") In [l + e'^^^")] . (6.11) 

Using standard methods p3|,^^ for the analysis of the TBA equations we find that their 
sum can be written as 



6 ^ r / e^"(^) \ 



(6.12) 



J x=—oo 
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where C{x) is Rogers dilogarithm function 



C(x] 



dt 



Int _^ln(l-t) 



l-t ' t 

Hence, the C25 + Ci is completely determined by the asymptotic behaviour of the solutions of 
(16.91) as X — > ±00: 



lim ipn{x) = In ((ri + 1)^ — l) 



lim ipn{x) 



In ((n- 25 + 1)2-1 

1 f sin^ (7rn/25+l) _ -, 
sin2{7r/25+l) ^ 



— 00 



for n > 2S 
ioi 1 < n < 2S 
for n = 1, 25* 



givmg 



C2S + Ci = 2- 



45 



(6.13) 



25+1 

independent of the doping (i.e. A). 

The individual values of the Cj are easily calculated for small and large doping corresponding 
to A and A ^ 00, respectively. In these cases the regions contributing to the integrals 
(|6.11|) are well separated and the functions v^„(x) take constant values in between. For small 
doping {A <ti 1) we find ipn{x) = for InA <^ vrx <^ with 



In {{n - 25 + 1)^ - 1) for n > 25 

< —00 for n = 2S 

l^(^j^piMs±m_i) forl<n<25 

sm"^(7r/2(S+l)) / — 



(6.14) 



Hence, near x ~ they behave as in the undoped system giving the central charge 35/ (5 + 1) of 
the 5f/(2)2s WZNW model. In the region around a; ~ InA the (pn<2S are solutions of the finite 
set of TEA of the minimal unitary model A4p [|2^ with central charge Ci = 1 — Q/{p{p + 1)) 



where p = 25 + 1 (this is the Ising model for 5 = 1 (see [TQ]), tricitical Ising model for 5 = 3/2, 
three-state Potts model for 5 = 2, tricritical three state Potts model for 5 = 5/2 and so forth). 
Putting everything together we find the leading contribution to the spin part ( |6.8| ) to the free 
energy at small doping 

ttT^ 35 7rT2 f 3 



fs 



1 



(6.15) 



Qv2s 5+1 6vi (5+l)(25+l) 
Proceeding in an analogeous way in the limit of large doping {A ^ 1) corresponding to a 
spin-(5 — 1/2) chain doped with spin-5 carriers we find 




for n > 1 
for n = 1 



(6.16) 
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for <^ TTx <^\nA. In this limit the low temperature contributions to fs can be written as the 
sum of a SU{2)i and a SU{2)2s~i WZNW model, the latter being the well known continuum 
limit of the pure spin-(S' — 1/2) Takhtajan-Babujian model: 

nT'QS-3 ttT^ _ ^^^^^^ 

For finite values of A the coefficients C2s and Ci in ( p.lOp have to be determined numerically. 
They are found to interpolate smoothly between their limiting values in ( |6.15| ) and ( |6.17|) . For 
S = 1 and S = 3 their doping dependence is shown in Figure 0. 



7 Summary and Conclusion 

To summarize, we have introduced a class of integrable models describing a magnetic system 
which upon doping interpolates between the integrable spin-S" and S — 1/2 Takhtajan-Babujian 
chains. These models arise when considering vertex models invariant under the action of the 
graded Lie algebra gl{2\l) with the local quantum spaces carrying the 'atypical' higher-spin 
representations [S]_^_. Their solution by means of the algrebraic Bethe Ansatz allows for a 
detailed study of their low temperature phase diagram. The spectrum of low lying excitations 
is described in terms of the dressed energies satisfying the TBA equations ([4.15|) and ([4.16|) . 



Without an external magnetic field the critical degrees of freedom separate into charge and 
magnetic modes as is well known in the Tomonaga-Luttinger liquid models for one- dimensional 



correlated electrons (see e.g. Refs. p8|,|29[|). Different from these models, however, one finds 
two branches of low lying modes in the magnetic sector which at small (large) doping can be 
identified with higher level SU{2)k WZW models and a minimal model (free boson). The WZW 
models have to be present in order to reproduce the well understood critical behaviour of the 
undoped and completely doped limiting cases. The second gapless magnetic mode, however, 
is quite peculiar: its appearence in the low energy of the undoped system is crucial to allow 
for the smooth crossover between the limiting cases ( 6.15| ) and ( |6.17|) subject to the constraint 
C2S + ci = const. 

The low-T behaviour of the S = 1 integrable model has motivated the proposition of an 
effective field theory of four (real) Majorana fermions as a possible starting point for studies of 



perturbations around the integrable model |10|. While free field representations could be used 



for the constituents of the undoped model, interaction terms between the two sectors had to be 
introduced to reproduce the change of the coefficients C2s and ci with the hole concentration 
observed in the exact solution. The possible form of this interaction term is constrained by the 
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SU (2)-symmetry of the model without a magnetic field. A similar construction of an SU (2)- 
invariant effective low energy field theory for the S > 1 models introduced here is possible 
by using the fact that the minimal models can be obtained within a GKO coset construction 



applied to |30, 31 



SU{2hs-i 8) SU{2), 

SUi2)2s ■ ^ • ^ 

In fact, the observed change in the conformal weights attributed to the magnetic modes be- 
tween the limiting cases of the undoped and the completely doped system appear to be just a 
'adiabatic' realization of this construction 

SU{2hs ® M2S+1 SU{2)2s-i ® 5f/(2)i . (7.2) 



one 



On the other hand, taking the limit H starting from the phase discussed in Section ^ 
may obtain a different field theoretical description of the S'f/(2)-symmetric phase: There the 
critical degrees of freedom can be described in terms of two free bosons each contributing c = 1 
to the sum C2s + Ci. For S = 1 this should give a complete description of the massless magnetic 
modes. It is likely that the apparent difference between the H ^ limit and the H = model 
can be understood as a rotation in the space of the effective fields (note that no physical field 
couples to one of the magnetic modes alone) [^. For S > 1 one has C25 + ci > 2 from (|6.13|) . 



Here, the difference to the finite field critical properties is similar to the one observed in the 
Takhtajan-Babujian models it is due to the appearence of gap for parafermionic degrees 
of freedom in the critical theory for any non-zero magnetic field. 
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A Two more Bethe Ansatze 



As mentioned above, the grading of the underlying algebra leads to different (though equivalent) 
sets of Bethe Ansatz equations (BAE) when using different orderings of the basis |TB[ . 

A second Bethe Ansatz can be obtained by starting from the pseudo vacuum \^s) by 
choosing a different reference state in the auxiliary eigenvalue problem for the amplitudes 
pa„-ai_ Following Refs. |13|,||,|16| the roots of 



1= II ^ T, a=l,...,Ni. 



are found to parametrize eigenstates of the transfer matrix ( p.l2|) with eigenvalues (Aj = Xj—iS, 
Va = — iS — i/2) 

jj, — Aj + t 



j=l ~ Aj 

a=l ^ I j=l /^-^j 



(The spectrum of the other transfer matrices defined in ( |3.5| ) can be obtained from this expres- 
sion by means of the fusion equations derived in Section ^). 

Alternatively, we may start from the fully polarized state |^^5_i ) = Y[n=i •5'— |, S'— 1)„ 

with maximal number of spin-^S* — |)-holes in the system. Now, eigenstates of (|2.12| ) obtained 
by adding particles to l^^^^i) and lowering the spin by are parametrized by solutions of 
the third set of BAE il 



« = 1, • • • , iVi • 

The corresponding eigenvalues are (Aj = Xj — iS — i, z/q, = z/q, — zS* — z/2 



A3 (^|{A,}f=i, m:^^) =-{^^ + 2^S + ^f H ^ ^ ' 



Xj- n 



+ + 2^Sf TT '"'^^^ + (/^ + 0" TT TT ^'""^^ • (A.4) 

a=i '^"-/^ i=i A,-/x t=i 
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For 5 = I Eqs. (|A.3|) become Lai's BAE for the supersymmetric t-J model [0,0. 



B Equivalence of the Bethe Ansatze 



In this appendix the equivalence of the sets (|2.17|) and (|A.3|) of Bethe Ansatz equations starting 



from the fully polarized state of spin-S* and S — 1/2 multiplets respectively is shown by means 
of a particle-hole transformation in the space of the rapidities. 

Following Refs. [^,|36| we rewrite the second set of Eqs. (|2.17|) as P^Ua) = with the 



polynomial 

k=l ^ k=l ^ 



According to the first of the + N^^ roots of P{uj) can be identified as Uo 

a = 1, . . . Nfi. Labelling the remaining A^^i ones as u'^ we have 



a: 



N 



-El^ f ^' 1^^! ) +ln(P(^.)/P(^p)) , (B.2) 



where Ca is a contour enclosing and z^^p = ± z/2 are the end points of the branch cut of 
the logarithm in ( |B.2| ). From the definition ( |B.1| ) we have 



I 

■) = T 

k=l 

which - when used in ( |B.2| ) - implies that 



P(A,±-) = T n i^i-^k±^) . (B.3) 



Using this identity in the first of Eqs. (|2 . 1 7| ) we obtain 



We continue by rewriting these equations as Q(A^) = with 



QiuJ) = iuJ + ^S)']JL-u'^-'-)-iuJ-^S)'YlL-u'^+'-) . (B.6) 

a=l ^ ^ 0=1 ^ ^ 
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Similar as above we can identify the first Nh + roots of this polynomial of degree L + 



with Xj and denote the remaining L — Nh = ones by A'. Proceeding as in 



we obtain 



Nh+Ni \ I i\ 

k=i V^a-Afc-a/ k=i 

fc=i 

Exponentiating this equation we obtain 



1 

27n 



dzln 



z + 



d^ 



InQ(z) 



7/' V * 



+ In ( Q{< + ^)/g(^:. - \ 



(.B.7) 



!3 



k=l ^fc 



K - 1 QK + 1) 



QK- 



fB.81 



Using this and the fact that z/^ solve the second set of Eqs. ( 2.17 ) together with the definition 
(|B.6|) the unprimed variables can be eliminated and we find 



n 



n 



Ai 



A'j, + ^ 



(B.9) 



Comparing Eqs. ( |B.5D and (B.9) with the Bethe Ansatz equations ( |A.3 ) the equivalence of 
the latter with ( p.l7| ) becomes evident. The proof of equivalence with ( |A.1| ) is completely 
analogous. 
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(a) 





X 



(b) 



Figure 1: (a) Schematic phase diagram of the doped spin chain in the ^-H plane: The bold 
line denotes the transition to a fully polarized state, interpolating between the saturation fields 
(|5.2|) for the spin S and S' = S*— 1/2 Takhtajan-Babujian chains. The left (right) shaded region 
corresponds to the undoped (completely doped) regime, (b) Phase diagram of the doped S = 1 
chain as a function of hole concentration x. 
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Figure 2: Dependency of Ci (dashed line) and C25 (full line) on the concentration x of holes for 
(a) 5 = 1, (b) 5" = 3. 
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